Compressed quantum simulation of the Ising model 
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In [R. Jozsa, B. Kraus, A. Miyake, J. Watrous, Proc. R. Soc. A 466, 809-830 (2010)] it has been 
shown that a match gate circuit running on n qubits can be compressed to a universal quantum 
computation on log(n) + 3 qubits. Here, we show how this compression can be employed to simulate 
the Ising interaction ol a ID-chain consisting out of n qubits using a universal quantum computer 
running on log(?i) qubits. We demonstrate how the adiabatic evolution can be realized on this 
exponentially smaller system and how the magnetization, which shows a quantum phase transition, 
can be measured. 
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The simulation of certain quantum system on a clas- 
sical computer seems to be an unfeasible task. The rea- 
son for this is the exponential growth of the required re- 
sources, like space and time, as a function of the number 
of considered quantum systems. As conjectured by Feyn- 
man [1] and proven by LLoyd [2], however, a quantum 
system can be used to simulate the behavior of another. 
The former one being such that the interactions between 
the systems are well-controllable and that the measure- 
ments can be performed sufficiently well. The suitabil- 
ity for the realization of such a quantum simulator has 
been shown for experimental schemes based on optical 
lattices or ion-traps [3-5] . Recently, several experiments 
using for instance trapped ions [6-9], neutral atoms [10] 
or NMR [11, 12] have realized quantum simulations. 

An important application of the quantum simulator is 
the study of the ground state properties of certain con- 
densed matter systems. Quantum spin models are well 
suited for the investigations of quantum phase transitions 
[13], which occur at zero temperature due to the change 
of some parameter, like the strength of the magnetic field, 
or pressure. The ID quantum Ising model, for instance, 
exhibits such a phase transition. It can be detected by 
measuring the magnetization as a function of the ratio 
between the interaction strength and the strength of the 
external magnetic field, which will be denoted by J here. 
Since the Ising interaction is relatively simple, this is a 
good model for the experimental demonstration of quan- 
tum simulation, even though it can be simulated classi- 
cally efficiently. 

One way to observe this quantum phase transition is to 
employ adiabatic evolution. The system, consisting out 
of n qubits, is initially prepared in the (simple) ground 
state of the Hamiltonian corresponding to J = 0. The 
adiabatic theorem tells us that if the parameter J is 
slowly increased the system will be at any time in the 
ground state of the Hamiltonian for this value of J [14]. 
Another way to transform the input state to the ground 
state of the quantum Ising model for the desired value 
of J is to apply a specific gate, U{J). In order to im- 
plement this gate, one usually approximates it using the 
Suzuki Trotter approximation [2]. Experimentally, the 
ground state properties, like the magnetization, M (J), 
of the Ising model for n = 2 [6] , n = 3 [7] and n = 9 [8] 
have been recently observed. 



Here we use a different approach, which makes use 
of the fact that certain quantum circuits, the so-called 
matchgate (MG) circuits (MGC), can be simulated by 
an exponentially smaller quantum system [15]. We ex- 
tend here this result and introduce new techniques to 
show that the evolution of the ID Ising model (including 
the measurement of the magnetization) of a spin chain 
consisting of n qubits can be simulated by a compressed 
algorithm running only on to = log(n) qubits. Here and 
in the following n is assumed to be a power of 2 and the 
logarithm is taken in base 2. More precisely, it is shown 
that the magnetization, M{ J), can be measured using the 
following algorithm: 1) First prepare the initial TO-qubit 
state p in = 1<8|ia) (yiU, with \ Vl ) = 1/V2(|0) - i |1)); 

2) evolve the system up to a certain value of J according 
to a specific unitary operator W(J) (see Eq. (7) below); 

3) measure the m— th qubit in y-direction, i.e. Y m . The 
expectation value of Y m coincides with the magnetiza- 
tion, M(J) (for n qubits) up to a factor — 1. The size of 
this circuit, i.e. the total number of single and two-qubit 
gates which are required to implement W(J), coincides 
with the one required to implement U(J) for the original 
circuit. Moreover, the error due to the Suzuki Trotter ap- 
proximation is the same as the one of the original system, 
since we are simulating the gates exactly. Due to the fact 
that this compressed quantum computation corresponds 
to the simulation of the Ising model not only the magne- 
tization, but also other quantities, like correlations can 
be measured. 

This result allows for the experimental measurement of 
the quantum phase transition of very large systems with 
current technology. Consider for instance, experiments 
with ion-traps or NMR quantum computing where say 
8 qubits can be well-controlled [8, 11, 16]. According to 
the results presented here such a system can be employed 
to simulate the interaction of 2 8 = 256 qubits. Of course, 
for such a large system the phase transition can be well 
observed. 

The outline of this paper is the following. First, we 
briefly recall the notion of MGs and some results related 
to MGCs. Then, we review some basic properties of the 
Ising model and its adiabatic simulation. After that we 
derive the MGs which are required for the adiabatic sim- 
ulation and show how the symmetry of the Ising interac- 
tion can be used to compress the whole simulation into 
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log(n) qubits. 

Throughout the paper we use the following notation. 
X, Y, Z denote the Pauli operators. If not stated differ- 
ently, the subscript of an operator will always denote the 
system it is acting on, or the system it is describing. Nor- 
malization factors as well as the tensor product symbol 
will be omitted whenever it does not cause any confusion 
and 1 will denote the identity operator. 

Let us begin by reviewing the notion of MGs. A MG 
is a two-qubit gate which is of the form A ® B, where 
A is acting on the even parity subspace, i.e. on span 
{|00) , 1 11)} and B is acting on the odd parity subspace, 
i.e. on span {|01),|10)}. Both, A and B are unitary 
and have the same determinant. A MGC is a quantum 
circuit which fulfills the following conditions: (i) the MGs 
act only on nearest neighbor (n.n.) qubits; (ii) the input 
state is any computational basis state (iii) the output is 
a final measurement in the computational basis on any 
single qubit. In [17] (see also [18, 19]) it has been shown 
that the output of any MGC can be classically efficiently 
simulated. In the following "n.n. MG" will be simply 
called "MG". 

In order to review how the efficient classical simula- 
tion is actually achieved we begin by introducing the 2n 
hermitian operators on n-qubits [20]: 

d=XI---I ••• c 2fc _i = Z-ZXI-I ••• 
c 2 =YI---I ••• c 2k =Z-ZYI-I 

Here, X and Y are in the fc-th slot for c 2 k-i and c 2fc , 
and k ranges from 1 to n. It can be easily verified that 
these operators satisfy the anti-commutation relations 
{cj,ci} = CjQ + cic 3 = 25jjl for j, I = 1, . . . , 2n. 

Let U be a MG acting on the qubits k and k + 1. Then 

2n 

tf Cj U = Y,Rjicu (i) 

2=1 

where R £ 50(2n, R) is a special orthogonal matrix with 
Rjl = 6 jt i VZ g {2k - 1, 2k, 2k + 1, 2k + 2} [18]. Consider 
now a MGC U = Un...U\ acting on the input state 
| 'I'm)- Let us denote by Ri the real orthogonal matrices 
associated to Ui (see Eq. (1)) and by R the orthogonal 
matrix associated to U, i.e. R = Rn ■ ■ ■ Ri- Noting that 
Zk = — «c 2 fc-ic 2 fc we obtain for the final Z measurement 
on qubit k 

(Z k ) = (* in | rfZ k U |* in ) - (2k - 1| RSR T \2k) , (2) 

where Sji = (* in | - iCjC t |* in ) for j ^ I and Sjj = [24]. 
For instance, for |Wj n ) = 10)®" the corresponding S ma- 
trix would be S 1 = i ©/j Yfe = 1 ® iY , where 1 denotes the 
(unnormalized) identity matrix on n dimensions. For any 
computational basis state S can be computed efficiently 
Hence, the simulation of the MGC, which consists then 
only of multiplications of 2n x 2n matrices, which can be 
efficiently determined, can be performed classically effi- 
ciently [25]. 



In [15] we showed that a MGC on n qubits can be com- 
pressed to a universal quantum computation running on 
log(n) + 3 qubits. If N denotes the size of the MGC, 
i.e. the number of single and two-qubit gates, then the 
size of the compressed algorithm is 0(N\og(n)). Thus, 
any polynomial-sized MGC, can be simulated by a uni- 
versal quantum computer of poly(n) size and exponen- 
tially compressed width O(logn). The idea was to apply 
the controlled gate, AU = |0) (0| <g> 1 + |1) (1| ® U, with 
U = S- 1 RSR T , to the input state |+) lo)® 106 ^. Mea- 
suring then the first system in X-direction leads to the 
desired result, (1| RSR T |2) [26]. Due to the fact that the 
required classical side-computation can be performed on 
log-space, the computation is indeed performed by the 
exponentially smaller quantum computer. 

Next, some basic properties of the Ising model are re- 
viewed. We consider the Hamiltonian, H(J) = Hq+JH\, 
to describe the ID Ising model with open boundary con- 
ditions. Here, H = Y^i=i %i an d H\ = Y^iZi Xi<E>X i+1 . 
In Appendix A the spectrum of the Hamiltonian and 
the corresponding non-degenerate ground state are re- 
derived [13, 20]. Since we want to compress the computa- 
tion we have to use here open boundary conditions, since, 
e.g. periodic boundary conditions would not correspond 
to n.n. MGs. 

The quantum Ising model shows a quantum phase 
transition at the critical value J — 1. At this point the 
second derivative of the magnetization is no longer con- 
tinuous. One way to experimentally observe this phase 
transition is to use adiabatic evolution [14] [27]. To de- 
scribe the evolution we discretize the Hamiltonian H(J) 
into L + 1 steps, as H(l), where I goes from to L and 
H(0) — H and H(L) = H + JmaxHi- In the case 
of the Ising Hamiltonian we have H(l) = Hq + J(l)Hi, 
with J(l) = J m axl/L. To obtain the ground state of the 
Hamiltonian H(J), the system has to evolve according 

to the unitary operator U(J) = YlfS'P Ui, where L(J) = 
JL/J max and Ui = e~ iH W At . Here, At = T/(L + 1), 
where T is the time of the evolution. The adiabatic limit 
is achieved for T, L — > oo (At — > 0). 

In a next step we use the Suzuki Trotter expansion 
to approximate Ui up to second order in At, i.e. Ui w 

e - t H At/2 e - l .J(l)AtH 1 At e - l H a At/2 + O ^ At 2y rj,^ wf) &p _ 

proximate the unitary U(J) by 

L(J) 

U(J) = II Ul ^o , (3) 
1=1 

where V = e - lH « At and U t = ViV with Vi = e - lJ ^ H ^ At . 

Hence, the procedure to simulate the adiabatic evo- 
lution and to observe the quantum phase transition is 
the following. First, the input state |0)®", which is 
the ground state of Hq, is prepared. Then the sys- 
tem is evolved according to the unitary U(J). After 
that, the magnetization as a function of J, M(J) = 
l/n£<*j|^kl*j}. where = U(J) |O) " denotes 

the ground state of the Hamiltonian H(J), is measured. 
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We will compress now the circuit explained above, 
which is running on n qubits to an equally efficient one 
which is running only on log(n) qubits. Note that for the 
standard classical simulation of a MGC, 2n dimensions 
are required, which corresponds to the Hilbert space of 
m = log(2n) = log(n) + 1 qubits (see Eq. (2)). How- 
ever, we will see that one can compress the algorithm 
even to log(n) qubits due to the symmetry of the Ising 
model. It should be noted here, that the compressed al- 
gorithm does indeed simulate the adiabatic evolution of 
the system. That is, it is not only possible to measure 
the magnetization using the compressed quantum com- 
putation, but also for instance correlation functions, like 
Xi ® Xi + \, could be measured. In this case, only the 
measurement would change, but the evolution would be 
the same. 

First we derive the rotations R i} which correspond 
to the MGs V = e- lHoAt and Vi = e - lJ ^ HlAt 
(see Eq. (1)). For a MG, U = e~ lHt we 
have H = i^^i^tCjCt, where h is antisym- 
metric. The corresponding rotation is then of the 
form R = e 4ht [18]. It is easy to verify that 
h = -l/2£Li(l 2fc - 1) (2k\ - \2k) (2k- 1|) and h x = 

-1/2 Efc^i(l 2fc ) ( 2fc + !| - |2fc + l)(2fc|). In order to 
make the connection to qubit systems we use the bi- 
nary notation, k = Ya=i 2 m_ % + 1 and write \k) — 
\k\, . . . , k m ). Then we have h n = i/21 ® Y m and there- 
fore R ee e 4h ° At = 1 ® e l2AtY ™ and R t ee e 4J ^ hlAt = 
(1 - cos(2J(/)At))(|l) (1| + \2n) (2n\) + cos(2J(/)Ai)l + 
sin(2J(l)At)J2kZl\2k + l) (2k\ - |2fc)(2fc + l|, with 
J (I) = Jmaxl/L. Thus, the compressed evolution 
which corresponds to the adiabatic evolution, U(J), 

is given by R(J) = y/RoTltl? T l^~\ with T i = 
RiRo and L(J) = JL/J max . Since the ground state 
of the Hamiltonian H(0), i.e. the initial state, is 
|0)®" the corresponding matrix S (see Eq. (2)), is 
given by S = il ® Y m . Putting everything to- 
gether we have (Z k ) = tv(U(J)^ Z k U(J) |0® n ) (0® n |) = 
{2k\ R{J)SR{J) T \2k - 1), where R(J) and S are defined 
above. Note that U(J) is acting on n qubits, whereas 
R(J) is only acting on m = log(n) + 1 qubits. In the 
following we will use a superscript to indicate the dimen- 
sion of the Hilbert space an operator is acting on. For 
instance the operators R™ , R™ and S m denote the oper- 
ators Ro,R t acting on U = ((F 2 )®" 1 , i.e. on the Hilbert 
space of m qubits. 

In order to determine the magnetization we can ei- 
ther measure all expectation values (Zk) and average 
over them, or, we use the fact that for large values 
of n the boundary conditions will not play an impor- 
tant role and that all expectation values will be ap- 
proximately equal to Z k ~ n / 2l like in the case of peri- 
odic boundary conditions. In this case, we would sim- 
ply measure only one expectation value. Here, we will 
compute the average since we are also dealing with sys- 
tems composed out of a few qubits only. To this end we 
use the fact that for any real antisymmetric matrix A 



we have Y2=\ ( 2k ~ l \ A \ 2k ) = l/2tr(Al ® iY m ) [28]. 
Since RSR has the properties of matrix A we have 
M(J) = i/(2n)ti(R(J)SR(J) T l <g> Y m ). Using now that 
i?o commutes with S = il® Y m and the orthogonality of 
the matrices we find 

M(J) = -l/(n)tr(R(J)p m R(J) T l ® Y m ), (4) 

with 

L(J) L(J) 

r(j) = n Ti = n r i r <> ( 5 ) 

i=i i=i 

and p in = l/2(S/i + 1) = 1 ® |yi) with \y x ) being 
the eigenstate to eigenvalue 1 of Y. 

In contrast to [15] we act here only on m, not m + 1 
qubits. The reason for this is that we do not apply the 
controlled unitary AS~ 1 R(J)SR(J) T , but we simply let 
the system evolve according to R(J). Note that this also 
allows for a simplification in the evolution and that this 
idea can be used for the simulation of any MGC. Due 
to Eq. (4), M(J) can be measured by preparing the m- 
qubit input state pt n = 1® |yi) (yi\ then applying R(J) 
for the desired value to J and measuring Y m . 

In Appendix B we show that due to the symmetry of 
the Ising model the algorithm can be even further com- 
pressed to rh = m — 1 = log(n) qubits. More precisely, it 
is shown that 

m(j) - ~tv(w(j) P tw(j)H® (6) 

where p"^ = l m -2<£> \yi) (yi\m anc ^ * ne ™ qubit operators 

L(J) 

w(j) = n °i( R t) T m 
i=i 

with Rf = 1 <g> e i2Aty ™, and Oi = U d (Rf l ) T '. Here, U d 
denotes the m-qubit phase gate which leaves all compu- 
tational basis states but \\)® m unchanged and the state 
|1)® A is mapped to e U2AT 

Thus, the whole quantum simulation can be performed 
on log(n) qubits. Explicitly, the compressed algorithm 
for the simulation of the Ising model reads: 1) prepare 
the initial m-qubit state pi n = 1 <g> \y\) 2) evolve 

the system up to a certain value of J according to the 
unitary operator W(J) [Eq. (7)]. 3) measure Y 7 %- The 
expectation value of Y^ equals —M(J) (for n qubits). In 
Fig. 1 (lower insert) the magnetization as a function of 
J for different values of the system size is shown. The 
phase transition can be clearly seen for values of rh > 5, 
i.e. for n > 32. To see how many steps in the Trotter 
expansion are required we compared in Fig. 1 (Fig. 2) 
the exact value of M(J) with the one obtained by the 
simulated evolution for rh = 3 (rh = 8). There, M(J) 
for different values of the total evolution time T, keeping 
At = 0.05 fixed, is shown. Considering the evolution for 
rh = 3 (Fig. 1) it can be seen that choosing T = 30, which 
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amounts to 600 steps in the Trotter expansion, suffices 
for the approximation. Recail that the same number of 
steps would be required for the simulation using the ex- 
ponentially larger system. However, the effect of errors 
will be much smaller in the compressed simulation. In 
order to demonstrate that, we consider in Fig 1 (upper 
insert) the error ( e IQZ )® lo s(") being applied after each 
step, OiRq in the compressed algorithm. The errors, a, 
are chosen randomly between and x = 10~ 2 (green line) 
and and x = 10~ 3 (red line). It should be noted that 
an error x < 10~ 3 basically does not affect the result. 




since we are simulating the gates exactly. However, in 
experiments, the error which will occur by implement- 
ing a gate on log(n) qubits instead of n qubits will af- 
fect the result less (see upper insert of Fig. 1 and Fig. 
2). Furthermore, the size of the circuit does not increase 
compared to the original adiabatic simulation. The (non- 




o i j 2 



FIG. 2: The same plot as in Fig. (1) but for rh = 8, i.e. 
N = 2 8 . The insert is plotted for T = 100, and L = 2000 
(At — 0.05), where the adiabatic evolution approximates the 
exact solution well. 



FIG. 1: The magnetization as a function of J for m = 3 
(i.e. n = 8). The solid line (online black) shows the exact 
value. The dashed-dotted line (red) shows the magnetization 
obtained for the total evolution time T = 10 and the number 
of Trotter steps, L — 200, (At = 0.05); the dashed line (green) 
for T = 30 and L = 600 (At = 0.05). The magnetization 
for T = 100 and L = 2000 (At = 0.05) basically coincides 
with the exact value. The upper insert shows M(J) for T = 
30, At = 0.05 (L = 600) solid line (blue) for different values 
of the error (see main text): lower line (green) for x — 10~ 2 
and dashed-dotted line (red) for x = 10 -3 . The lower insert 
shows the exact value of M(J) for rh — 1 (lowest line) to 
rh — 6 and for rh = 8 (most upper line). 



In summary, we have shown that the simulation of the 
Ising model of n qubits can be performed on a system 
consisting only out of log(n) qubits. The compressed 
quantum algorithm is an exact simulation of the adia- 
batic evolution: Whereas on the n-qubit system the uni- 
tary operation U(J) (see Eq. (3)) has to be applied for 
a certain value of J before the magnetization is mea- 
sured, the unitary operation W(J) (see Eq. (7)) has to 
be applied to the exponentially smaller system, before 
Yfh is measured. The error due to the Suzuki Trotter 
approximation is the same as for the original system, 



compressed) adiabatic evolution of the Ising model can be 
decomposed into the following MGs. First apply e lAtZi 
to each system i then the commuting two-qubit gates 
e ij(l)AtXi®x i+ i f Qr j Q ^ anc j afterwards the same gates 

for i even are applied. Thus, the number of required gates 
here is 0(n). Since W(J), is just a log(n)-qubit gate, it 
requires at most 0(n) basic operations. For instance the 
implementation of Rq = \®e~ l2AtYm , which corresponds 
to the n gates e lAtZi , requires only one gate. It is also 
interesting to compare the circuits of size m and m — 1 . 
One only needs to change R% to (i? " _1 ) T and to 
U ( i(R] n ~ 1 ) T ■ It is due to the phase gate, Ud, why the 
circuit cannot be further compressed. Finally, note that 
this compression of quantum simulation is not restricted 
to the Ising model, but can also be used for instance 
to simulate the XF-model, as already suggested in [15]. 
The results presented here can also be used to study and 
experimentally simulate quench induced quantum phase 
transitions, which would lead to new and interesting ob- 
servations [22]. 
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Appendix A: Exact diagonalization of H(J) 

Here, we briefly review a method to compute the eigen- 
values and the magnetization of the Ising model [13]. 
We introduce the fermionic operators a\. = Z...Za~^\ = 
l/2(c 2 fe-i-ic 2 fe) and a k = Z...Zo~\ = l/2(c 2k - 1 + ic 2k ). 
It is easy to see that those operators fulfill the anticom- 
mutation relations {ai, aj} = and {ai, a}-{ — Sijl. Since 

cr+ = (-^EUX^aJ we have Xi = (-l)£*=i a * afc (aj + 
ai) and Zi = 1 — 2a\ai. Thus, H = cftMa, where 
a = (ai, a\, . . . a n , a\) and 

/ A B \ 

B T A B 

M _ B T A B ... 

B T A B 

V B T A J 

with A = -Z and B = -J/2(Z + iY). Since 
(1„ <g> X)M(l n ® X) = —M the eigenvalues of M 
always occur in pairs, like Xi,— Xi- Diagonalizing 
this symmetric and real matrix M = UDU T , where 
D = diag(di, . . .d 2n ), with d\ = X n ,d 2 = -X n ,d 3 = 
X n -i,d4 = —X n -i,...d 2n -i = Xo,d 2n = — Xq, where 
X n > A„_i > . . . > Ao > and U an orthogonal matrix, 
which depends on J , we find H = WDb, with b = U^a. 
Hence, H(J) can be expressed in the new fermionic op- 
erators, b t as H(J) = Ya=i d 2 il+(d 2l -i - rf 2 ;)4;-i fo 2;-i 
with the ground state energy Y^i=i d>2i and the corre- 
sponding (non-degenerate) eigenstate, the vacuum states 
of the modes hi, i.e. = |0. . .0)^^. It is impor- 

tant to note here that the ground states energy is not 
degenerate. 

Let us now compute the magnetization M(J) = 
l/n(^{J)\Y^iZi |*(J))- Since Z t = -2ajai + 1 and 
ai = (ahi-i and a\ = (a) 2 i we have J2i a \ a i = 
a^Ra = tfWRUb, where R = £j \2l - 1) (21 - 1|. 

Thus, (*(J)\Ei4ai\*(J)) = (*(J)\$UiRUb\*(J)), 
and therefore M(J) = -2/nY^(?7l RU) 2l , 2l + 1, which 
shows a quantum phase transition at the critical point 
J = 1, as there the second derivative of M(J) is not 
continuous. 

Appendix B: Compressed circuit of width log(n) 

In order to prove that the Ising model can be simulated 
by a quantum circuit of width log(n) we introduce the 
following unitary operator: 

V , = l/v^(|0)<0| 1 ®l®Z m -»|l)<l|®l (8) 
i |0) (1| ® X®" 1 - 1 + i |1) (0| ® X® m ~ 2 ® Y m . 

This unitary can also be written as V = 
l/V2Etli a k\k)(k\ + p k \2n-k + l)(k\, where 
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a k = (-l) fe+1 and (3 k = (-l) fe Vfc < n and 
a k = Pk = —i Vfc > n. Below we show 
that Eg ee VRRV* = e -^tz 1 ®l®Y m = 
|0> <0j 1 <g> {R$) T + |1)(1| 1 <g> fl^- Furthermore, we 
have Ri ee VRiV* = |0) (0^ <g> O; + |1) (1^ <g> 0|, where 
= U d {Rf) T with E7 d = 1 + (e 42JAt - 1) |1) {if 
and O; = x^iOiYX®™. Note that both, 
i?o an d Ri are block diagonal, which implies 

that also R(J) = VR(J)V^ = UtHH^iRo) 
is block-diagonal. More precisely, we have 

R(j) = |o) <o| x ® njo i (i? ™) T ] + |1> (lli ® YliPiRf] = 

|0> (0^ <g> Wo + |l)(!|i ® Wj. It is straightforward 
to show that Wi = X®™VK *X®™. Using that 
VS^yt = — Zi (g> 5 m , we obtain for the magnetization 
M(J) = -l/(2n)tr((|0) (Od ® Wo^Wj + |1> <l| x <8> 

WiY^)^)) = -l/(n)tr((W l ® Y^W^l ® 
For the last equality we used the fact that 
tr(W l <8> FaWJi ® r A ) = tr(Wil ® Y^W^l ® 
which can be easily verified by noting that 
tr(W l ® ^mWjl ® F A ) is real [23]. This proves 
the validity of Eq. (6). 

The fact that Rq and R\ are block-diagonal can be 



easily seen as follows. We write i?o = cos(2At)l + 
sin(2At)E;l 2 1 |2fc)(2fc-l| 

|2n-2fc + l) (2n-2fc + 2| - h.c. Using that 

|2n-fc + l) = X® rn \k) it is easy to show that R = 
VR V^ = cos(2At)l - sin(2Ai) X)fe=i |2fc)(2fc- 1| - 
|2n-2fc+l)(2n-2fc + 2| - h.c. = e -OAtz 1 ®l®y m _ 
Similarly, we write Ri = cos(a;)ll + (1 — 

cos(x))(|l)(l| + |2n) (2n|)+sin(x) I^^i" 1 (|2fc + 1> <2fc| - 
|2n-2fc) (2n-2fc + l| - h.c.) + sin(ar)(|n + 1) (n\ - 
\n)(n + l\), where x = 2J(l)At. Then we 
find Ri = cos(x)l + (1 - cos(x)) |1) (1| + 
sin(^)Efe=r 1 (|2fc)(2fc+l| - |2fc+l)(2fc| + 

ism(x)\n)(n\ + X® m sin(ar) Efc^fl 2 *) (2k + 1| - 
|2fc + l) (2k\X® m . For any k such that 

fc < ti/2 — 1 |fc) = |0, k2, ■ ■ ■ , k m ) and therefore 
Ri = |0) (0| ® O; + |1) (1| ® O;', where ( = cos(x)l + 

(1 - cos(x))|l)(l| + sin(x)Efe=r 1 (|2fc)(2fc + l| - 
\2k + 1) (2k\ + ism(x) \n) (n\ and O' = x® m (0)*X® m . 
Note that Oi = Ud{R] n ) T ', where {7<j is the m-qubit 
phase gate, 1 + (e* - 1) |1) Thus, both, R and 

i?z are block-diagonal. 



